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Mott Transition vs Multicritical Phenomenon of Superconductivity and 
Antiferromagnetism - Application to k-(BEDT-TTF) 2 X - 
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Interplay between Mott transition and the multicritical phenomenon of d-wave superconductiv- 
ity (SC) and antiferromagnetism (AF) is studied theoretically. We describe the Mott transition, 
which is analogous to a liquid-gas phase transition, in terms of an Ising-type order parameter 
r/. We reveal possible mean- field phase diagrams produced by this interplay. Renormalization 
group analysis up to one-loop order gives flows of coupling constants, which in most cases lead 
to fluctuation-induced first-order phase transitions even when the SO (5) symmetry exists be- 
twen the SC and AF. Behaviors of various physical quantities around the Mott critical point are 
predicted. Experiments on ft-(BEDT-TTF)2X are discussed from this viewpoint. 

KEYWORDS: organic superconductivity, Mott transition, antiferromagnetism, renormalization group 



Emergence of magnetism and superconductivity (SC) 
around Mott transition 1 ) has been one of the central is- 
sues in strongly correlated electron systems. The Mott 
transition 1 ^ is referred to as a phase transition driven by 
a local Coulomb repulsion from metal to Mott insulator, 
i.e., an insulator with local spin moment, whether it pos- 
sesses an antiferromagnetic long-range order (AFLRO) 
or not. There are two distinct pictures for the Mott- 
insulating states: (I) In Slater's picture 2 ) justified in 
the weak correlation limit, the insulating state originates 
from doubling of a unit cell accompanied by the AFLRO. 
(II) In the Mott-Hubbard picture 1 ' 3 ) valid in the strong 
correlation limit, the system can be insulating without 
any symmetry breaking. This insulating feature arises 
from the Hubbard gap for the charge transfer due to a 
strong local Coulomb repulsion U. The superexchange 
interaction between localized spins causes an antiferro- 
magnetic phase transition at a temperature much lower 
than the Hubbard gap or U. Therefore, Mott transitions 
should be dealt with as a different degrees of freedom 
from antiferromagnetism (AF). 

The dynamical mean field theory (DMFT) has re- 
vealed that the Mott transition without any symmetry 
breaking occurs for the Hubbard model in infinite di- 
mensions. 4 ) A first-order Mott transition curve appears 
by varying a ratio of bandwidth W to U and terminates 
at a finite-temperature critical end point (Mott critical 
point). This Mott transition is analogous to the liquid- 
gas phase transition, as pointed out in a field theory. 5 ) 
Recently, by restoring spatial correlations neglected in 
the DMFT, the correlator projection method has given 
a similar structure of Mott transition in two dimensions 
but modified by spin correlations. 6 ) Diverging compress- 
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ibility is a common feature at the Mott critical point. 7 ' 8 ) 
To describe Mott transitions, one can introduce an 
Ising-type "order parameter" r] with its spatial variations 
as in the case of a lattice gas model. Here, we note that r] 
is not associated with any symmetry breaking. Actually, 
Kotliar et al. have developed Ginzburg-Landau (GL) 
theory of Mott transition in infinite dimensions where 
spatial correlations are absent. 9 ) For finite-temperature 
Mott transitions, the order parameter r] describes a static 
eigen mode of the dynamical mean field. Then, to discuss 
multicritical phnomena of the AF and the SC around the 
Mott transition, we need to introduce a scalar field 77, an 
0(2) vector field A for the SC and an 0(3) vector field 
m for the AF. 

Layered organic superconductor ^-(BEDT-TTF^X of- 
fers an ideal laboratory to address this issue. It exhibits 
a Mott transition as the effective ratio U/W is varied by 
changing X, isotope effect, and applying pressure. 10 ' 11 ) 
A phase diagram determined by recent experiments on 
the CI salt (X=Cu[N(CN) 2 ]Cl) is summarized in Fig. 1: 
There exists a first-order Mott transition curve separat- 
ing SC-AF phases at low temperature, which continues 
to higher temperature terminating at the Mott critical 
point (P cr ,T cr ) ^ (260 bar, 40 K). 11 15 ) The Mott crit- 
ical point has been observed with an anomaly in sound 
velocity as diverging charge compressibility ft e i 13 ^ which 
has been theoretically predicted. 4 ' T ' 8 ) This Mott transi- 
tion appears with a similarity to a liquid-gas phase tran- 
sition, as in other systems 1 ) and theoretical results. 4 6 ) 

In this compound, the Neel temperature Tn and 
the superconducting critical temperature T c meet at 
(Pbcn^bcr) — (280 bar, 13 K), showing the features of 
a bicritical phenomenon 10 12 ) (Fig. 1). GL theory for 
these two competing orders has been proposed, 16 ) and 
the analysis of the dynamical critical properties of NMR 
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Fig. 1. (Color) Experimental phase diagram of ^-(BEDT- 
TTF) 2 Cu[N(CN) 2 ]Cl with paramagnetic metal (PM), paramag- 
netic insulator (PI), antiferromagnetic insulator (AFI), inhomo- 
geneously coexistent region of AF and SC (SC/AF), and hystere- 
sis region of metal-insulator transition (MIT-HR). T c is deter- 
mined from AC susceptibility x /12 ^ and anomaly in ultrasonic ve- 
locity c, 13 ) while T N from NMR. 12 ) The first-order SC-AF phase 
boundary is determined from these three probes. 12 ' 13 ) Dotted 
curves represent spinodal instabilities for metastable AF (red) 
and SC (blue) determined from NMR, 12 ) x', 12 ) and anomaly in 
c. 13 ) The first-order MIT curve is determined from x' (green 
solid curve), 12 ) jump in resistivity p (black solid curve) 15 ) and 
maximum in dp/dP (open circles), 14 ) all of which lie in the 
dark shaded region. The MIT-HR is surrounded by spinodal 
instabilities (black broken curves) determined from hysteresis in 
p. 14 ) The Mott critical point is the end point (black filled circles 
for different experiments 12-15 )) of the observed first-order MIT 
curves. Crossover signaled by enhanced charge compressibility 
(green dashed curve) is determined from dip in c. 13 ) 
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1 revealed the SO (5) symmetry between the SC and 
which was first proposed for high-Tc cuprates. 17 ) 



The SO(5)-symmetric fixed point yields bicritical phe- 
nomena where order-disorder phase transitions become 
second-order. In contrast, away from the SO (5) sym- 
metry, a fluctuation-induced first-order phase transition 
from disordered to ordered phase 18 ) or a coexistent phase 
appears, depending on the strength of competition. 

Now, it is highly important to clarify the interplay 
between Mott transition and the competing two orders, 
i.e., AF and SC. In this paper, we develop a GL theory for 
this purpose. Global phase diagram is studied including 
the fluctuation up to one-loop order. Critical properties 
around the Mott critical point are also discussed. 

We start from the following free energy expansion in 
terms of three order parameters 77, A and m; 



r v r](r) 2 + (Vry(r)) 2 



2hrj(r) 



+r A A(r) 2 + (VA(r)) 2 
+r m rh(r) 2 + (Vm(r)) 2 
+u Am m(r) 2 A(r) 2 
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h plays a crucial role in describing the Mott critical point 
as introduced by Kotliar et al. in the original arguments 
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Fig. 2. Schematic MF phase diagrams. Bold and doubled lines 
represent second-order and first-order transitions, respectively, 
while circles the Mott critical point. 



of GL theory of Mott transition. 8 ) It is expected to obey 
h ~ ho(P — Pq(T)) around the Mott critical point where 
P controls (U/W) -1 and Pq(T) represents its value at 
the metal-insulator phase boundary, n represents a dis- 
tance from the bare critical point at temperature T^o 
(i = 77, A, m). It has the form, Ti ~ di(T — T^{P)). 
For SC and AF degrees of freedom, we take into account 
the effects of P on T^o by expanding it around the Mott 
critical point as T i0 (P) ~ T i0 (P (T v0 )) + h{P - P (T ? o)) 
with 6a > and b m < 0. Uij = 77, A, m) describes 
a quart ic coupling constant. g& (g m ) is a lowest-order 
coupling constant between 77 and A (m). 

First, we study the mean-field (MF) phase diagrams. 
If UAm ^> WAA,%m) the coexistence of SC and AF is 
prohibited. We also restrict ourselves to the case that at 
T < T cr , T c , Tn, the first-order Mott transition curve sep- 
arates the SC and the AF with Tn > T c , which is singnif- 
icantly important related to the ^-(BEDT-TTF^X sys- 
tem. Then, there are three possibilities as shown in 
Fig. 2: In (a), the first-order Mott transition runs out 
above Tc and Tn and ends at a critical point. If 
T cr ^> T c ,Tn, the average (77) takes (—)\/\^r j \/u r)r) in the 
metallic (insulating) side. Then, an effective value for 
T/\ — r m is modified by (g& + gm)(v)i which exhibits a 

2(#A + #m)vT^l7^ at tne 



jump given by 5rA-n 
Mott transition curve. In (b), the bicritical point of 
SC-AF transition coincides with the Mott critical point. 
Destiniy of this phase diagram is discussed below using 
a renormalization group (RG) analysis. In (c), the Mott 
critical point is inside the SC-AF phase boundary and 77 
is massive at the bicritical point. Then, the competition 
between AF and SC becomes the main issue. 16, 17 ) 

Next, we turn to the RG analysis of the fluctuations. 
We start with the case of Fig. 2(b), which is crucially 
modified by the fluctuation into Fig. 3 with the emer- 
gence of fluctuation-induced first-order phase transitions. 
The detailed discussion follows below. Introducing a 



Mott Transition v.s. Multicritical Phenomenon of Superconductivity and Antiferromagnetism - Application to tt-(BEDT-TTF)2X - 3 



Insulator 



Metal 



(a) 




(b) 



Fig. 3. RG phase diagrams modified by the fluctuations from the 
MF one given by Fig. 2 (b). 



mometum cutoff A = e A , gi = \JTdgiA £s and uij = 
I d u i:j A- £ with the integral I d = 2[(47r) d / 2 r(<i/2)]- 1 with 
£3 = 3 — | and e = 4 — d (d: the spatial dimension), we 
obtain the following RG equations for the dimensionless 
running coupling constants up to one-loop order; 

9 A = -^A + [(A r A + 2)^ A ^AA-^Vm^m^Am]-^i( 1 ) 

g'm -e3gm^[(N rn ^2)g m u rnrn -N A g A u Arn ]-gf r [2) 

(3) 



U'rjrj = -eUrjrj + + - (§% + 



u' AA = -su AA + [(N A + 8)u 2 AA +N m u 2 Arn ] +g\ (4) 
u' mrn = -eu mm + [(N m + 8)u 2 mm +N A u 2 Am \ +g m (5) 

^Am -^Am + f A 9m 

+ [(A r A+2)^ AA + (A^ m +2)^ mm +4^ Am ]^ Am , (6) 

where TVa = 2 (7V m = 3) is the dimension of A (fa) and 
' represents a derivative with respect to A. 

When #a = #m = 0, six fixed points were obtained in 
terms of e expansion. 19 ) In particular, there exist two 
nontrivial fixed points with nonzero iTs: 16 ) (A) O(N)- 
symmetric fixed point u AA = u m m = u Am = u* so with 
^f so = ivTs anc * N = ^A + A^m- (B) Biconical fixed 
point (wAA^mm^Am) ~ (0.09048, 0.08475, 0.05358) for 
d = 3. In our case of Aa = 2, N m = 3 and d = 3, 
the fixed point (A) is unstable except in a special case 
u AA Umm ~ while (B) is stable. 

With g Al g m > 0, however, there exist no nontrivial sta- 
ble fixed points, i.e., g A and g m are relevant couplings 
at the SO(5)-symmetric fixed point. Figure 4 shows 
the RG trajectories of running coupling constants with 
the SO(5)-symmteric initial condition u AA = Umm = 
UAm(= uq) with infinitestimal g A =g m - stay around 
uf so as long as g A ,g m ^ 0.5. (^a™ is n °t shown in Fig. 4, 
but it behaves similarly to u mm .) Therefore it is con- 
cluded that as long as g Al g m ^ 0.5, the SO(5)-scaling 
behavior is observed near the "bicritical point". When 
9a, 9m ^ 0.5, the instability of the SO(5)-symmetric 
fixed point is amplified and the bicritical scaling stops. 
u's flow to — oo, suggesting the fluctuation-induced first- 
order transitions. 18 ^ The difference among fi's rapidly 
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Fig. 4. RG trajectories around the SO(5)-symmetric fixed point 
obtained under the initial condition u^A — ^mm = UAm — 
at the scale g A =gm — ► 0. Curves in (b) correspond to different 
values of uq. 



increases in the run-away trajectories far from the SO (5)- 
symmetric fixed point (Figs. 4(a) and 4(b)), because of 
the difference between g's. Trajectories (g A ,g m ) shown 
in Fig. 4(c) depend on the initial value of uq = u AA — 
Umm — u Am . For the choices of uq > 0, g A becomes 
larger than g m due to the difference between N A = 2 
and N m = 3, and u AA becomes negative faster than 
Umm and u Am . Flows of are also similar to iiAA, 
indicating a first-order Mott transition. These RG re- 
sults suggest a phase diagram shown in Fig. 3(a) with 
the dominant SC over the AF. So far we have considered 
the RG flows around the SO (5) fixed point. On the other 
hand, if we consider the coupling constants far from the 
SO (5) fixed point and take g m ^> <?a, 9m remains to be 
larger than g A and u m m becomes negative faster than 
£tAA and u Am . Then, there occurs another possibility of 
phase diagram shown in Fig. 3(b) with the dominant AF 
over the SC. A phase diagram with a fluctuation-induced 
first-order phase transition to a coexistent phase is also 
possible. To determine which ordered phase appears, we 
must include higher-order couplings into the free energy 
and examine the parameters in detail. This is beyond 
the scope of this paper. 

We turn to the cases of MF phase diagrams, Fig. 2(a) 
and 2(c). In these cases, the r] field is massive around 
(Xbcn^Pbcr), i-e., | | / A§ ^ (Ao is the initial momen- 
tum cuttoff). (I) When I^I/Aq < 1, the fluctuation of 
r] is relevant and the RG analysis based on eqs. (l)-(6) 
applies for A^ = yj\r^\ < A < Ao- If an instability with 
negative fit's occurs before this scaling stops at A^, then 
the phase diagram is again given by Fig. 3(a) or 3(b). 
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Otherwize, the scaling for A < crosses over to that 
in the following case (II). (II) When \r v \/Al > 1, we can 
neglect the fluctuations of 77, and the problem is reduced 
to that of A and m only by putting 77 = ± ^ \ \ ju m in 
the case of Fig. 2(a) and 77 = in the case of Fig. 2(c). 
Furthermore, in the case of Fig. 2(a), the discontinu- 
ity SrA-m of ta — r m at the first-order Mott transition 
produces another scale for a crossover: The RG flows for 
uaa, ^mm, an d UAm with g's being zero in eqs. (4)-(6) 16 ) 
should be stopped at A c = ^5r&- m . Therefore, when 
5rA-m > Aq in Fig. 2(a), AF and SC are well separated 
by the first-order Mott transition, while the scaling the- 
ory in ref. 16 applies for the range of A c < A < Ao in the 
case of SrA-m < Aq. 

Now we discuss the phase diagram of ^-(BEDT- 
TTF) 2 X shown in Fig. 1. The Mott critical point is well 
above Tb cr , and the phase diagram looks like Fig. 2 (a). 
Therefore SC and AF occur independently separated by 
the first order Mott transition, and do not interact with 
each other. However Tn and T c almost meet (within 
1 K) at the "bicritical point", and the observed singu- 
lar behavior of NMR relaxation rate strongly suggests 
the SO (5) bicritical scaling. 16 ) These suggest that SC 
and AF communicate with each other across the Mott 
transition. This dilenma could be resolved by assuming 
small 5r A _ m /Ao, because the bicritical scaling governed 
by the SO (5) fixed point should be observed in a certain 
region before the coupling to rj grows. This small value 
of <fr* A _ m /Ao also leads to fluctuation-induced first-order 
phase transitions to SC and/or AF, as obtained in the 
RG analysis (Fig. 3). This is not excluded in Fig. 1, 
because the first-order Mott transition must be accom- 
panied by a hysteresis region which defines metastabil- 
ities of metallic and insulating states. The fluctuation- 
induced first-order phase transitions to SC and/or AF 
might be hidden in this hysteresis region, and further 
experiments on the metastability are required. More- 
over, if T cr approaches T^ CT in other materials by varying 
the geometrical frustration, 6 ) then the present scenario 
has a more chance to be observed experimentally. 

Lastly, we discuss critical properties around the Mott 
critical point (P cr ,T cr ). The critical properties are de- 
scribed by the Ising universality class. Then, the specific 
heat, the derivative of (77) with respect to h and its jump 
S(rj) obey the following scaling laws, 

Coc|T-T cr |- (7) 

d{n)/dh\ T , h=h cr oc (T - T cr )^ for T > T cr , (8) 

5(v)\t, h=h cr oc (T cr - Tf for T < T cr , (9) 

S(v)\ T =T cr ,h(x\h-h cr \^ s (10) 

with a ~ 0.013, 7 - 1.25, (3 ~ 0.31 and S ~ 5.0 in 
three dimensions. 21 ) Here, h cr is the critical value of h 
and h cr = if we take A = m = in our free energy. 
We note that the double occupancy D and U play roles 
of density and pressure, respectively, in the analogy to 
the liquid-gas phase transition. 4 6 ) Therefore, D and 
dD/dU exhibit the same critical properties as (77) and 
d(rj)/dh, and hence the local magnetic moment given by 
\ — y in the Hubbard model. The singularity at the crit- 



ical point and the discontinuity at the first-order phase 
boundary can also be observed in other physical quanti- 
ties such as resistivity, Drude weight, spin susceptibility, 
and longitudinal and transverse thermal conductivities. 
This is because these quantities can be expanded in terms 
of (77) /T and the first-order term exhibits a dominant sin- 
gularity. Then, in the vicinity of the Mott critical point, 
the same scaling properties as 77 can be observed by re- 
placing h with P. Especially, the thermal Hall conduc- 
tivity may be useful for the experimental scaling analysis, 
since phonon contributions are usually negligible. 

We have considered only the half- filling case, but one 
can consider the electron filling n as another controlling 
parameter in addition to U /W. In particular, the diver- 
gence of K e \ at the critical point 7 ' 8 ) has been observed as 
an anomaly of the sound velocity c in ultrasonic measure- 
ments by Fournier et a/., 13 ) although it is still difficult 
to estimate the scaling property from the experiment. 
The same anomaly for c can also be obtained with Ra- 
man scattering experiments. Extention of the present 
GL theory to include n is an interesting problem left for 
future studies. 
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